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The fixed point that governs the critical behavior of magnets described by the iV-vector chiral 
model under the physical values of N (N = 2, 3) is shown to be a stable focus both in two and three 
dimensions. Robust evidence in favor of this conclusion is obtained within the five-loop and six-loop 
renormalization-group analysis in fixed dimension. The spiral-like approach of the chiral fixed point 
results in unusual crossover and near-critical regimes that may imitate varying critical exponents 
seen in physical and computer experiments. 
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The critical behavior of the N- vector chiral model re- 
mains a controversial issue, being, at the same time, 
of prime interest for those studying phase transitions 
in helical magnets, stacked triangular antiferromagnets, 
and frustrated Josephson-junctions arrays both theo- 
retically and experimentally. In principle, this model 
should reproduce the main features of phase transitions 
into chiral states in above substances, but various ap- 
proaches and approximation schemes still yield quite 
different theoretical results (sec for review and 

guiiyjiKiipnn for most recent three_ and 

two-dimensional results). Even within the same, oth- 
erwise powerful field-theoretical renormalization-group 
(RG) approach the theoretical predictions are found to 
strongly, on the qualitative level, depend on the order of 
approximation and the reference point of the e-expansion 
The experimental data for the critical 
exponents obtained in several relevant materials turn out 
also to be considerably scattered (see p],^ppj p0j|?l| . f22| 
and references therein). 

In this work, we put forward the idea that the con- 
troversial situation one faces studying the phase tran- 
sitions in the aforementioned systems may reflect the 
quite unusual mode of critical behavior appropriate to 
the iV-vector chiral model under the physical values 
of N. To argue this conjecture, we analyze the criti- 
cal thermodynamics of the two-dimensional and three- 
dimensional chiral models within the five-loop and six- 
loop renormalization-group approximations in fixed di- 
mension and determine the structure of RG flows for 
N = 2,3. In the course of this study, the advanced 
resummation technique is used that proved to be effec- 
tive even near the borderline where the divergent RG 
series loose their Borel summability ||J2^,^J]. We will 
show that the critical behavior of the systems with planar 
and vector chiral ordering is governed by the stable fixed 
point which is a focus attracting RG trajectories in a 



spiral-like manner. Approaching the fixed point in a non- 
monotonic way, looking somewhat irregular, may result 
in a large variety of crossover and near-critical regimes 
that give rise to strongly scattered effective critical expo- 
nents seen in numerous experiments and computer sim- 
ulations. 

We start with the Landau-Ginzburg- Wilson Hamilto- 
nian of the D-dimcnsional chiral model 




where <p a , a = 1,2 are two iV-component vectors. Un- 
der vq > this Hamiltonian describes transitions into 
the phase with noncollinear ordering, i.e. into the chi- 
ral state. Since the lower-order renormalization-group 
approximations failed to give consistent results, we are 
in a position to study the critical behavior of the model 
(1) on the base of RG series of maximal length. For the 
three-dimensional chiral model such record, six-loop se- 
ries were calculated in ||, while in two dimensions only 
four-loop RG expansions were known JIo| . Aiming to 
obtain more reliable results, we extend the latter series 
to the five-loop order, using numerical values of the five- 
loop integrals evaluated in |25|| ; explicit expressions for 
the five-loop contributions to /3-functions and critical ex- 
ponents will be published elsewhere . 

The renormalization-group expansions are known to 
be divergent and some resummation procedure should 
be applied to extract proper physical predictions and nu- 
merical estimates. Taking into account the properties of 
Borel summability and the large order behavior of the 
series in two [ [To[ and three |^| dimensions, we resum the 
perturbative expressions by means of a Borel transforma- 
tion combined with an appropriate conformal mapping 
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|p7[ for the analytical extension of the Borel transform. 
With this resummation procedure we obtain the follow- 
ing approximants for each perturbative series R(u, v): 



E(R) p (a, b; u, v) = ^ B k (a, b; v/u) 



k=0 



x / dtt"e 



b„-t 



y(ut; v/u) 



[1 - y(ut;v/u)] 



- (2) 



where 



y(x;z) 



y/l- x/u b (z) - 1 
^\ —x/Tib(z) + 1 



(3) 



and Hb is the singularity of the Borel transform (evalu- 
ated in Ref . [p|Jl0| ) closest to the origin at fixed z = v/u. 
Varying two parameters a and b we are able to estimate 
the systematic error of the final results. We consider 
the 18 approximants for each (3 function that minimize 
the difference between perturbative orders: the approx- 
imants with b = 5, 7, . . . , 15 and a = 0, 1, 2 are used for 
the two-dimensional models [ fic| , whereas in three dimen- 
sions we choose the approximants with b = 3, 6 . . . , 18 
and a = 0,2,4 §. 

Applying this procedure to the evaluation of the com- 
mon zeros of the (5 functions, we find in three dimensions 
a chiral fixed point, as in In two dimensions the five- 
loop results also confirm the existence of the chiral fixed 
point found earlier at the four-loop level p0| . 

Then we look at the stability properties of these fixed 
points. Since the exponents w+ and w_ determining the 
fixed point stability are eigenvalues of the matrix 



<9Ar(u, v) d/3u(u, v) 



du 
d/3y(u, v) 
du 



dv 
dfiuju, v) 
dv 



(4) 



we consider the numerical derivatives of each pair of ap- 
proximants of the two /3- functions at their common zero. 
The majority of the eigenvalues, evaluated in two and 
three dimensions using 324 combinations of the approxi- 
mants for each case, has non-vanishing imaginary parts. 
Namely, for D = 2, N = 2 the complex eigenvalues are 
produced by 87% of working approximants, for D = 2, 
N = 3 - by 89%, and for D = 3, N = 2 - by 76.5%. 
For the physically important three-dimensional model 
with N = 3, all the approximants employed are found 
to give the eigenvalues with non-zero imaginary parts. 
To have reliable numerical estimates for u + and lo^ 1 we 
limit ourselves by the approximants that yield the chiral 
fixed point coordinates compatible with their final, prop- 
erly weighted values. The approximates leading to purely 
real eigenvalues of the stability matrix are discarded in 
the course of the averaging procedure. 

Thus, for the three-dimensional model with N = 2 we 
obtain 
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FIG. 1. RG flow for the physically important case D = 3, 
N = 2. Both /3-functions are represented by the approximants 
with a — 2 and 6 = 6. With this choice, the chiral fixed point 
coordinates (u,v) are (1.882,4.017). 



uj± = 1.40(35) ±i 1.00(45), 
u>± = 1.70(40) ±i 0.90(40), 



while for N 

u± = 
uj ± = 



0.85(20) ±i 0.95(15), 
1.00(20) ±i 0.80(25), 



5 loop, 

6 loop, 



5 loop, 

6 loop, 



In two dimensions we have for ./V 



lj± = 1.50(25) ±i 1.00(45), 4 loop, 
u ± = 2.05(35) ±i 0.80(55), 5 loop, 



and for N = 3 



u± = 1.30(25) ±i 0.50(35), 4 loop, 
u>± = 1.55(25) ±i 0.55(35), 5 loop. 



(5) 
(6) 



(7) 
(8) 



(9) 
(10) 



(11) 
(12) 



The appearance of imaginary parts comparable in mag- 
nitude with the real ones clearly indicates that the chiral 
critical behavior is controlled by the focus fixed point 
both in three and two dimensions |p8| . 

To further confirm this conclusion, we apply also the 
resummation by means of the Pade-Borel-Leroy tech- 
nique (see e.g. Refs. |3|]25||). The results are substantially 
equivalent, although the statistics is much less significant 
due to a large number of defective Pade approximants 
(i.e. possessing dangerous poles) generated during the 
resummation. 

In order to substantiate how the focus driven critical 
behavior manifests itself in the RG flows, we investigated 
the structure of these flows, given by the resummed five- 
loop (2D) and six-loop (3D) RG expansions. Here, we 
report four examples of RG flows corresponding to physi- 
cal values of N and generated, for certain values of D, by 
the typical working approximants. They are presented in 
Figs. 0, [|, ||, and ||, clearly demonstrating the spiral-like 
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0.5 1 1.5 2 2.5 3 



FIG. 2. RG flow for the physically interesting case D = 3, 
N = 3. Both /3-functions are given by approximants with 
a = 2 and 6 = 9, the chiral fixed point is located at 
(1.702,2.858). 
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FIG. 3. RG flow for the two-dimensional model with 
TV = 2. The approximants used are given by a = 1 and 
6 = 5 for fry. an d a = an< i 6 = 8 for In this case, the 
chiral fixed point is at (2.427, 5.045). 



approach of the chiral fixed point. Obviously, all the RG 
flows are quantitatively correct only within the regions 
where the singularity of the Borel transform closest to 
the origin is on the real negative axis (see fl|l0|); in Figs. 
0, H ||, and |J they correspond to unshaded areas. Nev- 
ertheless, we report also the flows in other parts of the 
coupling constant plane to present a complete qualitative 
picture @. 

From these figures another peculiar and interesting fea- 
ture emerges. For some manifold of the starting points, 
the RG trajectories have the coordinate v that grows very 
fast at the beginning and seem to reach the region of 
the first order phase transitions, but just before arriving 
there these trajectories drastically curve moving toward 
the stable chiral fixed point. Such RG flows may be con- 
sidered as a possible explanation of recent Monte Carlo 
simulations || in which a similar behavior was observed 
(in the space of Binder cumulants), but only up to the 
first rising up of the coupling v, interpreted as an evi- 
dence of fluctuation induced first order transition. 

Within the context of the results obtained, it is rea- 
sonable to discuss the fortune of the unstable anti-chiral 
fixed point. For almost all the working approximants, 
this point is seen in the region where the singularity of the 
Borel transform closest to the origin is on the real pos- 
itive axis, leaving the existence of the anti-chiral point 
doubtful. The fact that its location strongly oscillates 
with varying the approximants indicates that in this do- 
main of (u,v) plane the analysis is not robust. On the 
other hand, under the presence of the stable focus chiral 
point, the second, unstable fixed point is not topolog- 
ically needed on the separatrix dividing the regions of 
first and second order phase transitions. 

Finally we report in Fig. [s] the effective exponent v for 
N = 3 and d = 3, evaluated along the RG trajectories 
of Fig. 0. The exponent oscillates within a large range 



(0.47 < v < 0.72), even close to the stable fixed point. 
This provides a possible explanation of the scattered re- 
sults obtained in simulations and experiments; in fact 
the majority of numerical results, obtained when a con- 
tinuous phase transition was observed, are in this wide 
range. 

It is worthy to note that earlier the focus-like stable 
fixed points were found on the RG flow diagrams of the 
model describing critical behavior of liquid crystals |3C| l 
and of the 0(n)-symmetric systems undergoing first or- 
der phase transition close to the tricritical point |pl| . In 
these cases, however, the independent coupling constants 
had different scaling dimensionalities and played essen- 
tially different roles in forming the critical thermody- 
namics. Moreover, recently the focus driven chiral phase 
transition was observed in the three-dimensional model 
(1) within the three-loop approximation, but for quite 
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FIG. 4. RG flow for D = 2, N = 3. The approximants 
used are given by a = 1 and 6 = 5 for (3% and a = and 6 = 8 
for /3j-, the chiral fixed point coordinates are (2.201, 3.854). 
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FIG. 5. Effective exponent for iV = 3 and d = 3, along 
RG trajectories. Note that v is a multi-valued function of u, 
as a consequence of the non-monotonic behavior of the RG 
trajectories. 



unphysical values of N (N = 5, 6, 7) For real physi- 
cal systems having coupling constants of the same scal- 
ing dimensionality, the robust evidence in favor of phase 
transitions governed by the focus stable fixed point is 
presented for the first time. 
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